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Abstract 

Within the framework of QCD we compute renormalization constants for the 
strong coupling and the quark masses to four-loop order. We apply the DR scheme 
and put special emphasis on the additional couplings which have to be taken into 
account. This concerns the e-scalar-quark Yukawa coupling as well as the vertex 
containing four e-scalars. For a supersymmetric Yang Mills theory, we find, in 
contrast to a previous claim, that the evanescent Yukawa coupling equals the strong 
coupling constant through three loops as required by supersymmetry. 
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1 Introduction 

Dimensional regularisation [1,2], (DREG) is a remarkably elegant procedure which com- 
pletely dominates the radiative corrections industry associated with the standard model. 
Advocates of alternative regularisation methods rarely proceed beyond one loop (or ex- 
ceptionally two). The fundamental reason for the DREG hegemony is that (with little 
increase in caiculational difficulty) it preserves gauge invariance; that is to say, when the 
effective action is separated into a finite part (which is retained) and an "infinite" part 
(or more precisely, a part which tends to infinity in the limit that D = 4 — 2e — > 4) then 



the finite effective action satisfies the Ward identities of the gauge symmetry, without the 
necessity of introducing additional finite local counter-terms. 

DREG is, however, less well-suited for supersymmetric theories because invariance of a 
given action with respect to supersymmetric transformations only holds in general for 
specific values of the space-time dimension D. An elegant attempt to modify DREG so 
as to render it compatible with supersymmetry (SUSY) was made by Siegel [3]. The 
essential difference between Siegel's method (DREE0) and DREG is that the continuation 
from 4 to D dimensions is made by compactification, or dimensional reduction. Thus 
while the momentum (or space-time) integrals are D-dimensional in the usual way, the 
number of field components remains unchanged and consequently SUSY is undisturbed. 
(A pedagogical introduction to DRED was given by Capper et al [4].) 

As pointed out by Siegel himself, [5] there remain potential problems with DRED. One 
manifestation of this was demonstrated in Ref. [6], where it was shown that the variation 
5S of the action of a pure (no chiral matter) supersymmetric gauge theory is nonzero even 
with DRED. If SS gives a nonzero result when inserted in a Greens function this creates an 
apparent violation of supersymmetric Ward identities. With DREG this happens at one 
loop, but with DRED all explicit calculations to date have found zero for such insertions. 
For discussion see Refs. [7, 8] . 

We turn now to the application of DRED to non-supersymmetric theories. That DRED 
is a viable alternative to DREG in the non-supersymmetric case was claimed early on [4] . 
Subsequently it was adopted occasionally, motivated, for example, by the fact that Dirac 
matrix algebra is easier in four dimensions and in particular by the desire to use Fierz 
identities. One must, however, be very careful in applying DRED to non-supersymmetric 
theories because of the existence of evanescent couplings. These were first described [9] 
in 1979, and later independently by van Damme and 't Hooft [10]. After dimensional 
reduction the four-dimensional vector field can be decomposed into a D-dimensional vector 
field and a 2e-dimensional which transforms under gauge transformations as a scalar and 
is hence known as an e-scalar. Couplings involving the e-scalar are called evanescent 
couplings; in a non-supersymmetric theory they renormalise in a manner different from 
the "real" couplings with which we may be tempted to associate them. It has been 
conclusively demonstrated [11,12] that there exists a set of transformations whereby the 
/3-functions of a particular theory (calculated using DRED) may be related to the f3- 
functions of the same theory (calculated using DREG) by means of coupling constant 
reparametrisation. The evanescent couplings play a crucial role in this analysis, but in 
the literature on non-supersymmetric DRED their significance is often ignored, and there 
have been few calculations which explicitly take them into account. In a recent paper [13], 
four of us described in detail the evanescent coupling structure of QCD and calculated 
the gauge /^-function, {3, and the mass anomalous dimension 7 m through three loops. 
We found that at three loops j3 depends on the e-scalar Yukawa coupling g e , while 7 m 

1 Dimensional reduction in combination with modified minimal subtraction is often known as DR; we 
will use the terms DRED and DR interchangeably. 
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depends on both g e and the e-scalar quartic self couplings, A r . In this paper we extend 
these calculations to the four-loop level, when j3 also depends on A r . These results bring 
the precision of our knowledge of these quantities in DRED up to that in DREG [14-17]. 



2 Evanescent couplings 



The technical framework of our calculation is described in detail in Ref. [13]. Let us at 
this point emphasise once more the role of the evanescent couplings and in particular 
elaborate on the esse vertex. The part of the Lagrange density describing the latter is 
given by 

1 R 

c = •••- i E A ^ a6cd ^44 + --- , (i) 

r=l 

where the e a a are the e-scalar fields, and a, a' are 2e-dimensional indices. For the gauge 
group SU(N) the dimensionality R of the basis for rank-four tensors is given by R = 3 
for SU(2), R = 8 for SU(3) and R = 9 for SU(N), N > 4; for tensors Hf cd symmetric 
with respect to (a, b) and (c, d) exchange these numbers become R = 2, R = 3 and R — 4 
respectively [11]. We will restrict ourselves to SU(3); our basis choice reads 

j^abcd ^_ ^ jace jbde _|_ jade yfcce^ 

^jabcd ^ab^cd _|_ ^ac^bd _|_ ^ad^bc 

R abed = 1 (fiacfibd + flad^fcj _ ^ab^ed _ ^) 

Note that dimensional reduction of the original action yields Ai = g 2 , A 2 = A3 = 0, but, 
as we have already emphasised, this situation is not preserved by renormalisation. 

Once the tensors H® bcd are chosen, the Feynman rules are fixed in a unique way. It is 
straightforward to relate the result obtained from a different choice of the H^ bcd to each 
other. For example, for 

jjj-abcd ^^ab^cd 

2 

jjabcd = 1 ^ae § bd + gadgbcj ^ 

2 

jjabcd \_ ^jace jbde _|_ jade jbce'J 



one obtains 



Ai — A 3 , 
A2 = - ^A~i + 2A~ 2 

I ("Ai + A 2 ) . (4) 



A, 
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The renormalization constants for the evanescent couplings are defined through 

-0,a / ry £ a 
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where r g g e and r e£££ are the one-particle irreducible e-scalar-quark and four-e-scalar Green 
functions, respectively, the superscript "0" denotes bare quantities, and /x is the renor- 
malization scale. The charge renormalization constants are obtained from the following 
relations 
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with Z 2 being the wave function renormalization constant of the quark fields. 
Let us introduce the couplings 



n 9 ° 
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and define the corresponding (3 functions in the DR scheme: 



d/i Z 7T 



a 



DR 



a 



7T 
DR 



= — e- 



(3 e {a^ R ,a e ,{r] r }) = n 



7T 



+ 2— = , 

Z S DR I 5a 

/ ., r'ijklm 
i,j,k,l,m 

d a fl 



, a^dZ™ 
1 + 2^= — = , 

z dr 5a pR y 



1 



-1 



Q 



DR 



7T 



e- + 2-^ 

7T Z fi 



d/i 2 7T 

<9Z, 



1 + 2 



Z P <9a R 



e 

ijklm 



a 



DR 



i,j,k,l,m 



n 



Oe_y (r]iy (rj2\ l fm^ 71 

7T / V 7T 



7T 



7T 



7T 



d/i 2 7T 

d^A r ^dr <9Z Ar 



dr] r 



1 + 2 



Z\ r drj r 



i,j,k,l,m 



a 



DR 



7T 



(6) 



(7) 



(8) 



(9) 



(10) 



Here and in the following we do not explicitely display the dependence on the renormal- 
ization scale /i, i.e., a s = a s (fi) etc. Note that in the second line of Eq. ([8]), the 0(e) 
terms of /3 e and /3 Vr contribute to the finite part of 0^ R , and similarly for Eqs. ([9]) and 
( ITOl) . As we will see below, in order to compute the four-loop term of (3 DR one needs (3 e 
to two-loop and (3 Vr (r = 1, 2, 3) to one-loop order. 

In analogy to the (3 functions we introduce the quark mass anomalous dimension which 
is defined through 



da^ da e " d V „ 



i,j,k,l,m 



From this equation one can see that for the four-loop term of 7° R , the beta functions 
j3 e and (3 Vr are needed to three- and one-loop order, respectively, since the dependence of 
Z° R on a e (j] r ) starts at one-loop (three-loop) order [13]. 

The one-loop terms for /3 Vr and the three-loop expression for (3 e can be computed using 
standard techniques (see e.g. Ref. [18]), leading to the following non- vanishing coefficients 
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where rif is the number of active quark flavours and £3 = £(3) = 1.20206..., where ( is 
Riemann's zeta function. The one- and two- loop result of (3 e can be found in Ref. [13]. 



3 /3 e SYM and (3* YM to three loops 

In order to check our technical setup, we calculated the (3 functions of the strong and the 
evanescent coupling constant for a supersymmetric Yang Mills (SYM) theory. 

Such a theory is obtained from massless QCD when replacing the quarks by the SUSY 
partner of the gluon, the so-called gluino. Technically, this amounts to changing the 
colour-matrix for the fermion-gluon coupling of QCD from the fundamental to the adjoint 
representation of the gauge group. In addition, closed fermion loops have to be multiplied 
by an extra factor 1/2 in order to take into account the Majorana character of the gluino. 

SUSY requires that the gluino-gluon coupling a s equals the gluino-e-scalar coupling a e 
to all orders of perturbation theory, and therefore /3f YM = /3f YM for the respective (3 
functions. However, it was found in Ref. [19] that this equality is violated at the three- 
loop level. The result was interpreted [19,20] as a manifestation of explicit SUSY breaking 
by DRED, when employed in the component (as opposed to the superfield) formalism. 

Using the approach described above, we re-calculated /?f YM and /9f YM through three-loop 
order. When we set a e = a s , the result for /3f YM agrees with Ref. [21] (Ca is the Casimir 
of the adjoint representation of the gauge group): 



^SYM ' ( 1 - 

1 7T 



4 8 7T 64 V 7T 



+ 0(a 5 s ). (14) 



However, in contrast to Ref. [19], we find that 

P e = P s +0{a s ). (15) 

We draw the following conclusions from this result: (i) The expression quoted in Ref. [19] 
for the three-loop expression of j3 e is incorrect; considering the fact that this calculation has 
been done almost 25 years ago, it is probably impossible to trace the origin of the difference 
to our result; (ii) in a supersymmetric Yang Mills theory, the evanescent coupling constant 
a e renormalises in the same way as the strong coupling constant a s up to three-loop level, 
as required by SUSY; (iii) the setup of our calculation has passed a strong consistency 
check. 



4 /3 UK and 7^ K to four loops 

The direct way to compute the (3 S and 7 m function is based on the evaluation of the 
corresponding renormalization constants. For such a calculation one can exploit that the 
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divergent parts of a logarithmically divergent integral is independent of the masses and ex- 
ternal momenta. Thus the latter can be chosen in a convenient way (provided no infrared 
divergences are introduced). Up to three loops this procedure is quite well established 
and automated programs exist to perform such calculations (see, e.g., Refs. [22,23]). Also 
at four-loop order this approach is feasible, however, technically quite challenging. Thus 
we decided to adopt the indirect method discussed in Refs. [13,24]. It is based on the 
following formulae relating the quantities in DRED and DREG: 
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Let us in the following briefly discuss the order in perturbation theory up to which the 
individual building blocks are needed. Of course, the MS quantities are needed to four- 
loop order; they can be found in Refs. [14-17]. The dependence of a DR and m DR on a e 
starts at two- and one-loop order [13], respectively. Thus, (3 e is needed up to the three- 
loop level (cf. Eq. f[T2"|) ). On the other hand, both a DR and m DR depend on i] r starting 
from three loops and consequently only the one- loop term of /3 Vr enters in Eq. (Tl6l) . It is 
given in Eq. (TIB"]) . 

Two further new ingredients are needed for the four-loop analysis, namely, the three-loop 
relations between a DR and a^ s and between m DR and m MS . The two- loop results have 



already been presented in Ref. [13]. Parametrising the three- loop terms by 5a^ and 5m , 
we have 
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where the dots denote higher orders in a^ s , a e , and r] r . 5a^ and 5$ are obtained from 
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the finite parts of three-loop diagrams (see Ref. [13] for details). We find 
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We performed the corresponding calculation for arbitrary gauge parameter and use the 
independence of the MS-DR relation as a check of our result. Furthermore, let us stress 
that also the cancellation of the explicit In// 2 terms which occur at intermediate steps of 
the calculation is non-trivial. 
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Inserting Eq. (|T8|) into Eq. (TT6T) gives for the four-loop coefficients of the (3 function 
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where the four-loop MS coefficient is given in Eq. (8) of Ref. [14]. 
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Similarly we obtain for the four-loop coefficient of the mass anomalous dimension 
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where the four-loop MS coefficient 73 can be found in Eq. (6) of Ref. [15]. 



5 The four-loop supersymmetric case 

We saw in Section [3] that in the special case of SUSY, the relation a s = a e is preserved by 
the /3-functions through three loops. We now consider the supersymmetric case at the four 
loop level. The results in Section 0] were presented for the gauge group SU(3); however, 
we have evaluated those parts that are not related to the evanescent couplings 772 and 773 
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also for a general Lie group G. It is well-known that a simple change of color factors^] in 
addition to the statistical factor 1/2 for closed fermion loops, translates these terms into 
a supersymmetric Yang Mills theory. In this way, we can compare our four-loop results 
to the gauge /5-function (3f YM which was presented in 1998 [25]: 



0. 



SYM 
s 



TV 



3 n . 3 n 2 a s . 21 3 

4 8 n 64 



2 51 
f 128' 



7T 



+ 0{al 



[21) 



The method employed in Ref. [25] to obtain the four- loop result was very indirect, in 
particular relying on the NSVZ form [26, 27] of /3f YM . It is therefore a remarkable check 
on our calculations, and indeed those of Ref. [25], that we obtain precise agreement with 
Eq. (12T|) when we adapt our calculation to the supersymmetric described above. 

(As well as setting a e = a s we must of course set rjx = a s , ry 2 = % = 0). 

Turning to the mass anomalous dimension we have a similar powerful check. In the 
supersymmetric case the fermion mass term breaks SUSY; but j m (a.k.a. the gaugino 
/^-function) is in fact given in terms of /3 S by the simple equation [28]: 
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(This relationship^! between 7 m and (3 holds in both DRED and NSVZ schemes.) Through 
four loops we have at once from Eq. f)22p that 
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Quite remarkably, in the supersymmetric case we again find this agrees with our calcu- 
lations. This is again striking confirmation of our methodology and of the exact result 
Eq. (J22} • 



6 Conclusions 

In this paper we have applied DRED to QCD, and calculated both the gauge /^-function 
and the mass anomalous dimension to the four-loop level. These calculations required 
careful treatment of the evanescent Yukawa and quartic couplings of the e-scalar. In the 
supersymmetric limit we explicitly verified that the /3-function for the evanescent Yukawa 
coupling reproduces the gauge /3-function through three loops. 

2 Note that this procedure differs from the one outlined in Section [3l where we modified the Feynman 
rules in color space and re-evaluated the color factor for each diagram. Here, we simply replace Cf 
(fundamental Casimir) and T (fundamental trace normalization) by Ca (adjoint Casimir) in the final 
result. 

3 Note that our definition of 7 m in Eq. (fTTj) differs by a factor of two (and a factor of M) from the 
definition of (5m in Ref. [28]. 
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The popularity of the MSSM and the construction of the CERN Large Hardon Collider 
(LHC) has led to many increasingly precise calculations of sparticle production and decay 
processes, using DRED. The MSSM is a softly broken super symmetric theory, so we might 
well expect its dimensionless coupling sector to renormalise like the underlying super- 
symmetric theory, without worrying about evanescent couplings; to test this (in the same 
manner as described above) will require a generalisation of our calculations to incorporate 
scalar fields. In the MSSM, however, there is in fact one evanescent quantity which must 
certainly be considered: the e-scalar mass [29]. This exists also in QCD, but affects neither 
the gauge /3-function nor the fermion mass anomalous dimension on simple dimensional 
grounds so we have not considered it here. 

If, however, one wants to match MSSM calculations on to the Standard Model (or, for 
example, consider an intermediate energy theory produced by integrating out the squarks 
and sleptons [30]) then evidently the use of DRED inevitably means one must worry about 
evanescent couplings. Ref. [13] pointed out a couple of instances where naive treatment of 
the evanescent couplings has led to incorrect conclusions; we believe that careful treatment 
of the e-scalars in higher order calculations as presented in Ref. [13] and here will prove 
invaluable in matching MSSM calculations to low energy physics. 
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